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Abstract 

We investigate the relationship between the gap between the energy of the 
ground state and the first excited state and the decay of correlation functions in 
harmonic lattice systems. We prove that in gapped systems, the exponential decay 
of correlations follows for both the ground state and thermal states. Considering 
the converse direction, we show that an energy gap can follow from algebraic de- 
cay and always does for exponential decay. The underlying lattices are described 
as general graphs of not necessarily integer dimension, including translationally 
invariant instances of cubic lattices as special cases. Any local quadratic cou- 
plings in position and momentum coordinates are allowed for, leading to quasi-free 
(Gaussian) ground states. We make use of methods of deriving bounds to matrix 
functions of banded matrices corresponding to local interactions on general graphs. 
Finally, we give an explicit entanglement-area relationship in terms of the energy 
gap for arbitrary, not necessarily contiguous regions on lattices characterized by 
general graphs. 

1 Introduction 

On physical grounds, one expects that correlations in ground states of gapped many- 
body systems - quantum systems for which the energy of the lowest exitation is strictly 
larger than the ground state energy - decay exponentially. This implication of the 
energy gap to the exponential decay of two-point equal-time correlation functions is a 
well-established observation on non-critical quantum many-body systems ||T1I21|3]|31|5| 
|6||7]|8]|9). Yet, surprisingly perhaps, rigorous proofs of this implication for spin-models 
in higher dimensional cubic or general lattices have been found only very recently. 
Notably, Ref. |6) reexamines this question for finite-dimensional constituents, which 
has been generalized to the case of general lattices defined by graphs that may have 
any non-integer dimension in Refs. |3]|7). 

In this paper, we rigorously reconsider this question for a class of systems on gen- 
eral lattices the constituents of which are infinite-dimensional quantum systems: the 
quasi-free case of harmonic systems on general lattices (see Figure^. This harmonic 
case is particularly transparent, as the entire discussion of state properties can be done 
in terms of the second moments of the states. This case has notably taken a central 
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Figure 1 : A general graph with a distinguished region I <Z L (gray area, dark gray 
oscillators). Lines indicate the edge set, the number of bold lines is the surface area 
s(I). 



role in recent discussions of questions of entanglement scaling in many-body systems, 
as a class of natural physical systems where sophisticated questions on the scaling of 
entanglement - in form of, e.g., entanglement-area relationships 1 1 01 1 1 lllT2l and other 
issues 11311141 [T51I16II17I - are relatively accessible. 

The general lattices are not assumed to be necessarily translationally invariant, and 
the underlying graphs may have any spatial dimension. Harmonic systems, needless to 
say, have been considered many times before, in particular the case of nearest-neighbor 
interaction 1181 : Here, we investigate in a very comprehensive manner the implications 
of a gap to the two-point correlation functions and the converse direction for all har- 
monic lattice systems with local Hamiltonians reflecting general couplings in position 
as well as in momentum, both for the ground state and Gibbs states (non-zero tem- 
perature). Notably, these findings enable us to formulate an entanglement-area law for 
regions of arbitrary shape within the considered context of harmonic systems on gen- 
eral graphs. Such harmonic systems model discrete versions of Klein Gordon fields, 
vibrational modes of crystal lattices or of ions in a trap, or serve as approximations to 
not strictly harmonic systems. In contrast to, for example Ref. @, we will not rely on 
using variants of Lieb-Robinson 1 19 1 bounds on the group velocity to prove the validity 
of our bounds. Instead we will make use of generalizations of the methods introduced 
in Ref. 1201 to the case of general lattices (see Figure|2j. 

• We briefly introduce the notion of quasi-free states as ground or thermal states 
of harmonic Hamiltonians on generic lattices described by graphs, with coupling 
in position and momentum coordinates. 

• For local harmonic Hamiltonians on generic lattices we prove that whenever the 
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Exponential decay 




Figure 2: Schematic summary of obtained results. 



system is gapped, exponential decay of the correlations in the ground state with 
the canonical coordinates follows (Theorem 1). 

• Conversely, for systems (not necessarily locally) coupled in position, we demon- 
strate that for sufficiently fast decay of the correlations, the existence of a spectral 
gap can be deduced. This includes sufficiently strongly algebraically decaying 
correlations (Theorem 2,3). 

• These findings give rise to a proven equivalence of exponentially decaying cou- 
plings and the existence of a spectral gap for systems coupled in position on 
generic lattices (Corollary 1). 

• We prove that for gapped harmonic lattice systems on generic lattices in Gibbs 
states (thermal states), exponential decay of the correlations follows (Theorem 
4). 

• We give a connection 122111 lllT2l to entanglement properties in many -body sys- 
tems, in this instance of harmonic lattice systems. We show that similarly to 
Ref. HI 21 . the entropy of a subsystem of a gapped lattice system is bounded from 
above by a quantity linear in the boundary area of the distinguished subsystem, 
which may be of arbitrary shape (Theorem 5). This generalizes the area theorem 
of Refs. 11 1 111 21 to systems on general lattices. 

• We finally discuss several examples in detail, such as instances of disordered 
systems or thermal states of rotating wave Hamiltonians. We also investigate a 
case of non-local algebraically decaying interactions. 

Note that simultaneously, similar findings were published in Ref. 1211 . where in con- 
trast to this paper, emphasis was put on ground state properties of translationally in- 
variant harmonic systems on cubic lattices, but where in turn, non-local interactions 
and critical cases were considered in much more detail. 
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2 Considered models and main results 



2.1 Harmonic systems on general lattices 

We consider quadratic bosonic Hamiltonians on generic lattices. This means that the 
lattice is characterized by general graphs [23 1 . The set L of vertices with cardinality | L \ 
is associated with the set of sites in the lattice, each of which corresponding to a bosonic 
degree of freedom. The simple graph G = (L, E) is identified with the lattice, where 
E is the edge set or, equivalently, the adjacency matrix characterizing a neighborhood 
relation between physical systems. In a cubic lattice, say, G will represent just this 
general lattice, regardless of the interaction, which may be finite-ranged beyond nearest 
or next-to-nearest neighbor interactions. 

We will consider paths in the sense of standard graph theory, as sequence of vertices 
connecting a start and an end vertex, each of which are connected by an edge. For two 
vertices i, j G L, the integer dist(i, j) denotes the graph-theoretical distance, so the 
length of the shortest path connecting i and j 1231 . with respect to the edge set E. 

(i) The Hamiltonian is assumed to be of the form 

H = fV p p + x T V x x. (1) 

where V X ,V P G ]Rl L l x l L l are positive matrices, with x T — (aci, and 
p T = (pi, ...,P\l\) being the canonical coordinates satisfying the canonical com- 
mutation relations. Hence, we allow for any general coupling in position and 
momentum coordinates (but no coupling involving both simultaneously). We 
write (V x , p )ij, i,j = 1, \L\, to label the element of the matrices V x , v that 
belongs to the coupling between two vertices i, j, respectively. For simplicity of 
notation, these couplings in position and momentum can be subsumed into the 
tuple C = (G, V x , V p ). Whenever two sites i, j G L are coupled in position, the 
element (V x )ij will be non-zero, and similarly for the momentum coordinates. 
The range of the interaction is taken with respect to the metric dist(-, •), such that 
V x and V p inherit a neighborhood relation. 

Note that V x and V p can in turn be conceived as adjacency matrices of a weighted graph 
with the same vertex set L, but a different edge set. The above description is similar 
to the assessment of the correlation function for generic spin systems in Refs. Q]|7]- 
Note that such bosonic harmonic lattice systems resemble to some extend the concept 
of graph states for spin or qubit systems in the sense of Refs. 1241 1251 . compare also 
Refs. tT2lfl4lfT6l. 

The dimension of the underlying lattice G can take any positive value in the fol- 
lowing sense: We may define a sphere S r (i) for some i € L, centered at site i with 
integer radius r as 

S r {i) :={l e L : dist(7, i)=r}. (2) 

Then there exists a smallest d > of the lattice, notably not necessarily integer fTl . 
such that for all r G N, 

Bup|5 r (<)| < cr d - 1 (3) 
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for some c > 0. This number d is taken as the dimension of the lattice. Note that 
for cubic lattices this dimension coincides with the natural underlying spatial dimen- 
sion, and the above expression is the same for open and periodic boundary conditions. 
Similarly we can define a ball centered at site i £ L with radius r as 

B r (i) :={l 6L: dist(/, i) < r} (4) 

and can express the maximal volume 

r 

u d , r := sup |B r (0| ^cYy- 1 (5) 

of a ball with radius r in the graph theoretical sense. For a consideration of the di- 
mension of general graphs and the discussion of self-similarity in this context, see, 
e.g., Ref. 1261 . For our purposes, the dimension will only enter the bounds through the 
volume of a ball of some radius. 

Concerning the locality, we assume the following: 

(ii) Most Hamiltonians that we consider are local, corresponding to a finite-ranged 
interaction, if not otherwise specified. This means that there exists an m G N, 
such that for alH G L 

(V x )i,j = 0, (VpJu = (6) 

for all j G L for which 

dist(i, j) > m/2. (7) 

Note that this in turn means that the Hamiltonian is of the form as in Eq. and 
can be written as 

XcL, 
diam(X)<m/2 

X C L being subsets the diameter of which satisfy 

diampf) = sup dist(i, j) < m/2. (9) 



It will turn out to be convenient to collect the \L\ conjugate pairs of canonical 
coordinates in a vector f T = (x\, ah^i ,£>i, ...,p\n), the entries of which satisfying 
the canonical commutation relations (CCR), giving rise to a symplectic scalar product 
1271 . The Hamiltonian 



H = fT { o v P J f ' (10) 

can now be brought into diagonal form by means of a symplectic transformation, so by 
means of linear transformations S G Sp(2\L\, R) preserving the symplectic form 1271 . 
defined by the skew-symmetric symplectic matrix 



V x 



1 

-1 



(11) 
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Let O G be the orthogonal matrix that brings Vx V p V x l into diagonal form, 

O t V x 1/2 V p V x 1/2 =:D = diag(rfx, • • ■ , d\ L \). Then 

( v 1/2 o \ 
s=l v ' Q ° d2) 



and with f = Sf \ the Hamiltonian finally takes the form 

I 

(2d] /2 y/* 



F = 2^(^ + 1/2), h:=^0^. (B) 



»=i 

The ground state energy is thus given by 

Eo = f>! /2 = tr[(vV*V p VVr /2 ) = tr[(V x V p )y% (14) 
i=i 

The energy of the first excited state is given by 

Eh. = E + AE, (15) 
A-E := 2A^ 2 (FVV p ^ )=2A V2 (T4 y p)! (16) 

where we denote by A m i n (yl) the minimum eigenvalue of a matrix A. The correspond- 
ing ground state of such a Hamiltonian is a quasi-free (Gaussian) state 1271 . completely 
characterized by the first and second moments. Quasi-free bosonic or Gaussian state 
means that the characteristic function - expectation value of the Weyl displacement 
operator - is a Gaussian function in state space. Equivalently, a state is Gaussian if 
its Wigner function is a Gaussian function in state space. The second moments can be 
embodied in the covariance matrix with respect to the ground state 



{W + W- (1 7 ) 



= ({(r' i -{f! i )),(r' j -(m+) 

w + m 

Here, first moments vanish as the ground state is the vacuum. We find 

y = D 1/2 ffir 1/2 , (is) 

and in the original coordinates we have 

7 = SiS T (19) 

= (v-v\v^v p vv*?i*v-v*) © (yw<ywv r vw)-wvw) . 

This expression simplifies significantly in case of commuting V x and V p . Often, V p = 
1, in which case 

7 = V- 1 ' 2 © (20) 

see Ref. 1101 . The language chosen here is the one used also in the assessment of 
entanglement properties of harmonic chains 11011131 or more general harmonic systems 
on lattices II llll2ll27lll4l[T5ll28l . mildly generalized to coupling in both position and 
momentum and to general lattices. 
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2.2 Exponential decay of correlations in ground states 



The first result links the spectral gap of the Hamiltonian to the exponentially decaying 
correlation functions of the system. Note that the bound on the right hand side depends 
only on the dimension of the underlying graph, the gap, the parameter characterizing 
the range of the interactions, and bounds to the coupling strength. 

Theorem 1 (Exponentially decaying correlation functions) Consider a Hamiltonian 
on a general lattice of dimension d with a coupling C = (G, V x , V p ) of finite range m 
as defined above. Then for all i, j € L with dist(i, j) > m the ground state satisfies 

\(xiXj)\ < K\\V p \\exp[-dist(i,j)/{\, (21) 
|<M')| < K\\V x \\exp[-dist(i,j)/Z], (22) 

where 

_ \\V x V P \\ 1/2 v d , m /2 _ 2m 

1 /2 

where || • || denotes the operator norm andAE = 2\^{V X V P ) is the energy difference 
between the first excited state and the ground state. 

Note that if the energy is expressed in units of ||V^V^||, then the correlation length 
depends only on the gap and the range of the interactions. 



2.3 Spectral gap from algebraically decaying correlation functions 

In this subsection, we consider an instance of the converse direction, compare Ref. 
l2l : we assume that the system has algebraically decaying correlation functions in the 
position coordinates and assume that we have a - not necessarily local - coupling in 
these coordinates and V p = 1, then we can conclude that the system must be gapped if 
the decay is sufficiently strong. This is made more rigorous in the subsequent theorem. 

Theorem 2 (Spectral gap from algebraic decay) Consider a sequence of couplings 
C {n) = ( G W ] y(n) } n £ ~N, on general lattices = (L^, fiW) of dimension 
d^ n \ Let K > 0, Kq > 0, 77 > supd*-™) =: d, and c := supc*-™-* < 00 as defined in 
Eq. (fJJ. If the ground states satisfy 

\(4 n) ^ n) )\ < [ k j, ,. .... M ' '■ (24) 



\Kdist(i,j)-" fori + 3, 



for all i,j G L'™' and for all n, then the energy difference between the first excited and 
the ground state satisfies 

inf AE^ =: AE > 7 > 0, (25) 

" K Q + cK((l + 77 - d) 

where Q is the Riemann zeta function. 
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Note that, if in addition to the assumptions of the above theorem the coupling V^ n > 
is (i) of finite range m for all n, and (iia) the momentum correlations decay as in Eqs. 
( 1241 or alternatively, (iib) the coupling has a finite coupling strength, sup ||Vw|| < 
oo, 1 it immediately follows from Theorem 1 that in fact the correlations decay ex- 
ponentially. This is an interesting observation, showing that in a sense the class of 
harmonic lattice systems is not rich enough to - roughly speaking - show all "types of 
decays" of the correlations. 



2.4 Equivalence of spectral gap and exponentially decaying corre- 
lation functions 

For exponentially decaying correlations we can always conclude that the system is 
gapped. This statement can be formulated as follows. 

Theorem 3 (Spectral gap from exponential decay) Let K > and £ > 0, and con- 
sider a sequence of couplings C^ n ' = (G^ n ' , V^ n \ 1), n € N, on general lattices 
= (iW^W) of dimension S- n \ d := swpd^ < oo, and c := supc (n) < oo. 
If the ground state satisfies 

\(4 n) xP)\<Kexp[-dist(i,j)/Z) (26) 

for all L^ n ' and all n, then the energy difference between the first excited and the 

ground state satisfies 

inf AE^ = :AE> — - > 0, (27) 

A(l + cLti- d (c Vs)) 

where Li\-d is the polylogarithm of degree 1 — d. 

Together with Theorem 1 this establishes for locally coupled systems the following 
equivalence. 

Corollary 1 (Equivalence of spectral gap and exponentially decaying correlations) 

Consider the sequence of couplings C'™* 1 = (G'™' , V W , 1), n € N, of finite range m 
on general lattices G^ n > = {L^ n \ of dimension dS n > with d := supii'™' < oo, 

and c := supc^ < oo. Then the following statements are equivalent. 

( i) There exist constants K,£ > such that both correlations with respect to the 
ground state satisfy 

\{xi n h^)\,\(p\ n) p^)\ < Kcxp [-dist(i,j)/$ (28) 

for all i,j £ L^ n ' and all n. 

( ii) The energy difference between the first excited state and the ground state satisfies 

inf A_E ( ") =: AE > (29) 
and the coupling is of finite strength, sup |j V( n ) j < oo. 



Note that (iib) follows from (iia) using the same methods as in the proof for Theorem 2. 
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This result establishes that in local harmonic systems on generic lattices, being non- 
critical in the sense of being gapped is equivalent with finding exponentially decaying 
correlation functions. 



2.5 Expontial decay in Gibbs states 

In this subsection, we consider thermal Gibbs states 1271 corresponding to some tem- 



perature T > 0, 



Q(T) = 



exp(-#/T) 



(30) 



tr[exp(-7J/T)] 

Such states are again quasi-free (Gaussian) states and can thus also be uniquely char- 
acterized by their covariance matrices, appropriately modified for this case of non-zero 
temperature. For a study of two-point correlations in Fermi-systems at non-zero tem- 
perature see also Ref. In the diagonal basis dl3t . it takes the form 



7(T)' 



£1/2 



(31) 



1 + 2(exp(2 J D 1/2 /T) - 1) 1 

© D~ 1/2 \l + 2(exp(2D 1/2 /7 1 ) - 

and in original coordinates we have 

7 (T) = S 7 (T)'S T = 7 (0) + V-^ G (V^V P V^) l/ V" 1 / 2 

®V^{V^V P V^Y 1/2 GV^A32) 



G 



:= 2 (exp (2(V^V p V^y/T) - 1 



(33) 



In the following we will now assume that matrices V x and V p commute. This yields a 
simplified covariance matrix 



7CO 
G 



7(0) + (V~^X /2g ) ® (V x 1/2 V p ~ 1/2 G), 



(cxp(2(14^ P ) 1/2 /T) 



1 



(34) 
(35) 



Furthermore, we require the following assumption on the lattice. This assumption is 
very similar to the ones in Ref. 0j, and is satisfied for a large class of natural lattice 
systems. 

Assumption 1 (Lattice structure) Consider a Hamiltonian on a general lattice G = 
(L, E) with couplings C = (G, V x , V p ) as in Theorem 1. Then, it is assumed that there 
exist constants Iq>Q and v > such that 



exp(—^idist(i, fc)) exp(—fxdist(k, j)) < Iq exp(— vdist{i, j)), 

keL 

for all i, j G L and for 

/i := - 



(36) 



V !c V p \\-(AE/2)- 



m 



(37) 
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For example, for a cubic lattice in d dimensions, we have that L = [1, n] xd , and 
dist(i, j) = Yls=i ~ 3s\' Therefore, we arrive at 



n 



-/xdist(i,fe) _— /Ltdist(fc,j) 



IIE 



e — *«l e — j« I 



n 

5=1 



-/i|i s -fc s | 



i& -js\ + 



e 2 ^ - 1 



5=1 fe 5 =l 
e -A*(i5+i5-2) i g— M(2n— is— is) 



< 



2 

//c 



3 2 M _ j 



e 2 M _ 1 



exp(- M dist(i,j)/2), (38) 



i.e., in this cubic case, we have v = p/2. 

Under the previous assumption, one can arrive at the subsequent statement on two- 
point correlations in Gibbs states. Note that it is not merely trivially true that the 
correlations are shorter-ranged in thermal states: there are examples of Hamiltonians 
where Gibbs states at higher temperatures have longer-ranged correlations than at zero 
temperature, see Example|2] 

Theorem 4 (Correlations in systems at finite temperature) Consider a finite-ranged 
Hamiltonian corresponding to C = (G, V x ,V p ), [V x ,V p ] = 0, on a general lattice 
G = (L, E) with assumptions as in Theorem^\and equipped with Assumption 1. Then, 
for dist(i, j) > m, the Gibbs state with respect to some temperature T satisfies 



where 



K(T) := 



\(xiXj)\ < K(T)\\V P \\ cxp [-dist{i,j)/£\ , 
|<Pift>| < K(T)\\V x \\exp[-dist(i,j)/Z], 



\\V x V p \\^v d 



m/2 



{AE/2f 



4l \\V x V p \\/(AE/2f 



exp 



AE / I (A£/2)2 

4P4^ P || 



1/2 



- 1 



(39) 
(40) 



(41) 



max \ — , 

" log 



2m 



VxV P \\ 

\\-(AE/2) 2 



(42) 



This theorem ends the list of main statements on the relation between the system 
being gapped and the decay of the two-point correlation functions with respect to the 
ground state and thermal states. In the next section we will study the implications of the 
above findings on the entanglement scaling of distinguished regions of general lattices 



2.6 Entanglement scaling in ground states of harmonic systems on 
general lattices 

The above statements on the decay of correlation functions have immediate implica- 
tions on the scaling of entanglement. More specifically, in a lattice system we may 
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distinguish a certain part I C L and ask for the degree of entanglement between the 
degrees of freedom of this region and the rest of the lattice. This question is that of 
the scaling of the geometric entropy of this region. This question goes back to semi- 
nal numerical investigations in Refs. 113 II 1321 (see also Ref. 1331 ). suggesting a linear 
relationship between the geometric entropy and the boundary area of a distinguished 
region. In turn, in Refs. II 1 1 II 21 a rigorous relationship between the boundary area and 
the entropy of the region in cubic harmonic lattice systems has first been analytically 
established. 

In fact, in the light of the above findings, one may infer the validity of such an 'area 
theorem' for general lattices, not only for cubic lattices, exploiting exactly the same 
methods of proof. Defining the surface area of a distinguished region / of the whole 
lattice as 

s ( J ) : = E E 1 ' (43) 

i£L\I jel 

dist(ij) = l 

see Figure^ one arrives at a bound to the von-Neumann entropy, 

E I s = S{Qi) = -tt[g J logM}, (44) 

of the reduced ground state, gj = tr^y [g], with respect to the distinguished region. 

Theorem 5 (Entanglement-area law for general lattices) Consider a general lattice 
G = (L, E) of dimension d equipped with a coupling (G, V, 1). If the coupling is of 
finite range m, the entropy of entanglement satisfies 

Es ~ log(2)(A£/2)2 S(7) ' (45) 

where Li\-2d is the polylogarithm of degree 1 — 2d, AE the energy gap above the 
ground state, s(I) the surface area of I, c as defined in Eq. (fJJ, and 

t:= , "1 ^ (46) 



log 



[V 



In this form, the contribution of a geometrical factor, as well as one originating from 
the correlation length, is very transparent. Note that, using Theorem 4 and the bounds 
derived in Ref. 1121 . one arrives at a similar bound for the distillable entanglement of 
thermal states. 

Indeed, this theorem covers the entanglement-area relationship in all generality 
for gapped harmonic models on generic lattices. In particular, for cubic regions with 
volume L d in d dimensions this means that the geometric entropy is bounded by ex- 
pressions linear in L d_1 . For critical fermionic quasi-free systems, in turn, there are 
instances where one finds a behavior of L d_1 log L for the geometric entropy 1341 1351 . 
In contrast, there is numerical evidence that for bosonic harmonic systems, even in the 
critical case, the validity of the entanglement-area relationship of L d_1 is preserved 
1361 1371 . To strictly prove or refute this relationship in this critical case constitutes one 
of the intriguing open problems in the field. 
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3 Proofs 



This section will contain the proofs of the statements made before. We will make 
extended use of the handedness of the interaction matrix in the sense of the metric 
dist(- , ■). The main ingredient will be polynomial approximations to matrix functions 
of matrices. We generalize a statement of Ref. |2()| to the case of general lattices. This 
extends the generalization put forth in Ref. 1121 for general cubic lattices. Note that for 
statements of the type following, random walk methods would similarly be suitable. In 
particular, for the case of nearest-neighbor interaction, this analysis has been done, see, 
e.g., Ref. j 18 1 for a comprehensive introduction. 

In the following we will be needing the subsequent lemma which is concerned with 
the range of matrices in the sense of dist(-, ■). Note that the matrix power is taken in 
the sense of the ordinary matrix power. 

Lemma 1 (Range of matrix powers) Let the \L\ x \L\ matrix Abe a coupling matrix 
on a general lattice G = (L, E). If A is of finite range m, i.e., Ai_j = Ofor dist(i,j) > 
m/2, m S N, then A n , n£N, is of range nm. 

Proof. This can be proven by induction over n. For n = 1 we have A n = A is of 
range nm. Now suppose A n has range nm, i.e., we have {A n )ij = for dist(i, j) > 
nm/2. Now, 



keL 

Suppose dist(i, j) > m(n + l)/2. It follows that for dist(i, k) > nm/2 we have 
{A n )i,k = as A n is of range nm. For dist(z, k) < nm/2 we have 

m(n+ l)/2 < dist(i, j) < dist(i, fe) + dist(fc, j) < nm/2+ dist(fc, j), (48) 

i.e., dist(fc, j) > m/2, i.e., Af. j — 0. In turn, for dist(fc, j) < m/2 we have 

m(n+ l)/2 < dist(i, j) < dist(i, k) + dist(fc,j) < dist(i, k) + m/2, (49) 

i.e., dist(i, k) > nm/2, i.e., (A n ) i j { = as A n is of range nm. Finally, we hence 
arrive at (A^ n+1 ^)i t j = for dist(i, j) > m(n + l)/2, which concludes the proof. □ 




(47) 



3.1 Proof of Theorem 1 



From the covariance matrix 7, Eq. dl9> . we have 




(50) 



(51) 



12 



We can now compute the matrix (V X ^ 2 V P V X ^ 2 ) l//2 and similarly (V P ^ 2 V X V P ^ 2 ) 
using the power-series expansion of the square root, which is valid for |ar| < 1, 



(1-x) 



(52) 



k=l 



where the coefficients are bounded by < 1 

-1/2 



(yt"v p vy*y = \\v x v p \\- l/2 (i - (i - v,y 2 v p v^ 2 /\\v x v p \\)) 

/ oo 

= IIWII" 1/2 [i + J2a k o(i-D/\\v w v p \\) k o 



-1/2 



k=l 



Now, 



i.e., 



D = 



X 1/2 v p v x ^o = (y-^oy 1 v p v x [v-^o 



(vy%v^y 1/2 = \\v x v p \\-v 2 vv 2 w p , x v-v 2 , 



W p , x := \1 + Y, a k^-{V P V X )/\\V X V P 



k=l 



and similarly 



(vy*v x vy*) 1/2 = \\v x v p \\-^x/ 2 w x , p v p ^ 2 . 



(53) 

(54) 

(55) 
(56) 

(57) 



Note that clearly, {Vx^VpV^ 2 ) -1 ^ 2 is a positive matrix, whereas W p . x is in general 
not a symmetric matrix (and the matrix power is the fc-fold concatenation of matrix 
multiplication). Thus, the correlation functions take the form 



= \\v x v p \\-^ 2 (v p w p , x ) itj , 
(p^) = \\v x v p \\-v 2 (V^^ . 



(58) 
(59) 



Assuming V x and V p to be of range m/2, i.e., (V x )ij = (V p )ij = for dist(i, j) > 
m/2, we have that V X V P and V P V X are of range to. Note that for a matrix to be of finite 
range in the sense of the above Lemma it is not required that it is symmetric. Hence 
we can conclude that (V x V p ) k and (V p V x ) k are of range km. Thus 



(W P , X ) 



hi+ E *k((i-(y J y*)/\\v x v p \\) k ) t 

fc>[dist(i,j)/m] 



(60) 
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For the diagonal terms we have \(W p , x ) hl \ < \\W P , X \\ = \\V x V p \\^ 2 /(AE/2) and for 
the off-diagonal terms we find 

(W P , x ) id \ < a k \\l~(V p V x )/\\V x V p \\\\ k 

fc> [dist(ij") / m] 

<- " t (>-W)'- 

k>im(i,j)/m-\ v 11 xVpl1 / 
The same bound holds for the entries of W X , P , i.e., for all i,j € L we have 

< p^y a _ (AwV is,(y)/m . (62) 



(W Pi:c ). . , (W XtP ) 



(ae/2) 2 y \\V X V P \\ 

For the correlation function {x%Xj) this yields 

...... < ll*^ll 1/a y,, TM (AE/2)Y St(M/TO 



dist(fc,z)<m/2 

Now dist(i, j) < dist(i,fc) + dist(fc, j) < to/2 + dist(fc,j), and thus we find for 
dist(i, j) > to 

irMI < \\V X V P \\ 1/2 sr inn \(^ KV' St(,J,/W 



feSL, 
dist(fe,i)<m/2 



< UK, 



keL, 
dist(fc,i)<m/2 



For | (piPj) | one can argue in the same manner. Noting that the sum on the right hand 
side is just \B m /2(i)\, we finally have for dist(i, j) > to 

|( ^- >| * (aw w rg' ' (65) 



\\V x V p \\ 1/2 v d , m /2 llT , „ A (Ai?/2) 2 \ dis,(,j)/(2m) 



'<**>' * w '^'^wJ ' (66) 

which concludes the proof of Theorem 1 . □ 

3.2 Proof of Theorem 2 

For each n, the energy gap is given by 

AEW = 2A min ((y(")) 1 / 2 ) = ||( ^ ( „ ) 2 ) _ 1/2|| . (67) 
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The operator norm is bounded from above my the norm corresponding to the maximum 
sum of elements of any row of the matrix, and hence we may write 

||(^(«))-l/2|| < sup J- | (( y(n))-l/2 kj | (68) 

< sup [Kq + K dist(i,j)- r > 

and hence 

||(T/(«))-i/2|| < sup ( Ko + KY^r^ J2 M' ( 69 > 

ieL< " ) V r=l dist(t,j)=r / 

and therefore 

|| (y(»))-Va|| < sup (tf + #fVW)IJ 

= ifo+c (T,) ifC(l + r)-rf W ), (70) 

where £ is the Riemann zeta function, and c- n > > is defined as in Eq. l|3}. Thus, we 
can conclude that 

/ \ 2 2 

~ K + d n )K((l +Tj- " Kb + cK"C(l +»? - <0 ' 
independent of n. □ 



3.3 Proof of Theorem 3 

For each n, the energy gap is given by Eq. J67> . Analogous to the previous proof, we 
have that 

|| (1/ (»))-i/2|| < K sup [i + yW-r/OI&WI ] (72) 

< A ^1 + c< n) £ exp(-r/0r d( " ) - 1 J 

< A(l + cLii-aic- 1 ^)) < co, 

where is the polylogarithm of degree 1 — d. As the above bound holds for all n, 

we can conclude that 

inf A£ (n) > — 2 , rr— > 0. (73) 

" K(l + cLii_ d (e-V£)) 

□ 
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3.4 Proof of Theorem 4 



In the finite-temperature case, the co variance matrix 7(T) is given by 

7 (T) = 7 (0) + (Vp(W)- 1/a G) © (k(K^)" 1/2 g) ■ (74) 

We will now proceed by demonstrating the exponential decay of the entries of the 
matrices G and (V X V P )-^ 2 = (V X /2 V P V X /2 )-^ 2 (V x and V p are assumed to be com- 
muting) and will then apply the assumption on the lattice to finally arrive at Theorem 
2. To show the exponential decay, we will now proof a generalization of a theorem of 
Ref. 1201 to matrices reflecting general lattices. This latter work is concerned with the 
exponential decay of entries of matrix functions of banded matrices. The proof will 
be very similar to the one in Ref. |20|, only that the notion of a distance is different: 
the distance from the main diagonal of a matrix versus dist(-, •) in the graph theoretical 
sense. So quite surprisingly, the ideas of Ref. |20| carry over to this more general case 
with little modifications. Note that a mere naive embedding of the potential matrices 
into banded matrices would be insufficient to find the above conclusion, as for every 
ordering of physical systems, the range of the banded matrices could not be kept finite. 
We nevertheless state the full proof here for completeness. 
We denote by 

a:=A min (n b:=\\V\\ (75) 

the minimal and maximal eigenvalue of a real symmetric matrix V of finite range. Let a 
function / : (D — > (D be such that fotjj is analytic in the interior of an ellipse 2 e x , x > 1> 
with focii in —1 and 1 and continuous on e x . Furthermore suppose (/ o ijj)(z) G R for 
z G R. Here 

^:€->C ^)= (6 -° )g +0 + 6 - (^) 



We can now state the generalization of Ref. |20|. 

Theorem 6 (Exponential decay of entries of matrix functions) Let V = (Vij) be a 

positive real symmetric matrix of finite range m, Vij = for dist(i,j) > m/2, let 
f : (D — > (D be such that it fulfills the above assumptions. Then there exist constants K 
and q, < K, < q < 1 such that 

\[f(V)] h3 \<Kq dis '^\ (77) 

where 

K := max(||/(y)||,-^max|(/oVO(z)|), (78) 



2/ 1 



<?:=(-] (79) 



2 We denote its half axes by a and fl, a > 1, /3 > 0, a > /3. It is then completely specified by the 
parameter % = a + j3. Note that, if a is known so is (3 as 1 = a 2 — 1 . 
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Proof. For a function g : (D — > C analytic in the interior of the ellipse e x , % > 1, 
continuous on e x , and with g(z) <G R for real z, one has 

infi max \g(z) - p(z)\ : p e P k \ < 2 max\g(z)\, (80) 

[-1<z<1 J x fc( x _l) zG£x 

where denotes the set of all polynomials with real coefficients and degree less than 
or equal to k. This result is due to Bernstein (see Ref. 1291 for a proof) and applies to 
the function / o ip. As V is assumed to be of range m, we have that (V)) is of 

range km for all polynomials pefl.. Thus, for dist(i, j) > km/2 we have 

(KV _1 (^))kj=o. (8i) 

For given i ^ j we now choose k such that 

k = [2dist(i,i)/m] - 1, (82) 
i.e., we have k < 2dist(i, j)/m < k + 1, which yields 

\(f(V)k 3 \ = \(tf°1>)('>l>- 1 (V)))i*-(py>- 1 (V)))ij\ (83) 

< ik/o^)^- 1 ^))-^- 1 ^))!! 

max \(f o ip)(z) - p(z)\ , 

—1<Z<1 

where the last equation follows from the fact that the spectrum of is contained 

in the interval [—1,1]. Applying Bernstein's theorem, we know that there exists a 
sequence of polynomials p(„) of degree k that satisfy 

lim max |(/ o i/j)(z) - pr n \(z)\ 

n— *oo — 1<z<1 

= inf | max J(/ o tp)(z) - : p E Pk 

2 

- ~( max | (/ o | 

X {x - 1) zee x 

„ / i \ 2dist(ij)/m 

<—X-max\(foil>)(z)\ - . (84) 

where the last inequality follows from the choice of k and x > 1. Along with the fact 
that (/OOk-s < 11/(^)11 this concludes the proof. □ 



We will now apply Theorem 3 to the function / defined as 

= (85) 

- reflecting the matrix function that is needed in order to evaluate the second moments 
of a Gibbs state - and the matrix V X V P in order to prove Theorem 2. In the notation of 
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Theorem 3 we have V = V X V P , which is of range 2m, a = (AE/2) 2 , and b = ||14VJ,| 
The function f(z) is analytic for all 



z G <C\{z G R : z < 0} 



(86) 



and therefore / o ip is analytic for all z G C with 0?(z) < (a + b)/(b — a). So for 
X = b/(b — a) = a + f3 the function / fulfills the requirements of Theorem 3 and 
|(/ o tp)\ attains its maximum at z = —a. Thus, 



2.V 



max |(/ o ijj)(z)\ > 



X - 1 *6 

To summarize, we have 

4||F x F p ||/(A£/2) 2 



> 



e 2(a(l-a/(46)))i/2 /T _ j 

2 



(87) 



ll/WII 



ex pl n 



(AB/2) 



1/2 



1 - 



- 1 



(AE/2) 2 

\\v x v p \\ 



dist(i,j) / m 



(88) 



To bound the entries of (V^VJ,) -1 / 2 we could also apply Theorem 3, however, for this 
special function, the bound can be given in a more straight-forward manner. To this 
end, consider again the power series expansion of the square root, Eq. ( 1521 . Denote by 
a and b the minimal and maximal eigenvalue of V — V X V P , respectively. As matrices 
V x and V p are assumed to be of range m, V is of range 2m. Thus, for given i ^ j 
choosing k = [dist(z, j)/m~\ — 1, i.e., k < dist(i, j)/m < k + 1, yields 



( v - m ), 



< 



- (^t + Y^ar(l-V/b) k ^j 

iv/b)- 1 ' 2 - + j2 mi - y i h ) k ^j 



i 

Vb 
1 

Vb 
1 

—= max 

y/b a/b<x<l 



a r(l-a;) 

r=k+l 



oo 

<-j= J2 *r(l-a/b) k 



r=k+l 



Vb 



a\fc+i Vb 



— !-=■) <— !-=■ 



a s dist(i,j)/m 



(89) 



For diagonal terms, we have | (V 1 ^ 2 )i,i\ < \\V 1 / 2 || = I /{AE/2), i.e., for all i, j we 
have 

\2\ dist(i,j)/m 



((^^ P )- 1/2 



'■;/ 



< 



(AE/2) 2 



1 



(AE/2)' 
\\V X V P 



and thus, together with the assumption on the lattice, 

U \\V x V p \\ 3 / 2 /(AE/2) 4 



{{V x V p )- l ^G) iti 



< 



-dist(i ,j)v 



cxp AE ( 1 



(AE/2)' 

4||y x V p | 



1/2 



/T 



(90) 



(91) 
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The last step is now analogous to the proof in the zero temperature case, Eq. (1641 . We 
finally arrive at 

< ^\\V x V P \\ 3/ W m/ 2\\V P \\/(AE/2f c _ Mi)/2 (92) 
" ex P (A£(l-^) 1/2 /T)-l 

< 4/ ||y^ p f/^ >m/2 |i^||/(A^/2) 4 c „, dlsta?)/2 (93) 



{y p (v x v p )-v*G 
(v x (v x v p )-v 2 g) 



cxp 

for dist(i, j) > m, which concludes the proof. □ 

3.5 Proof of Theorem 5 

From entanglement theory we know that an upper bound to the entropy of entangle- 
ment is given by the logarithmic negativity 1 38 39 40l l41ll42ll43l . which in turn can 
be bounded from above by an expression depending only on the two-point position 
correlation function with respect to the ground state 1 12 1, 

j£L\I 



From the proof of Theorem 4, Eq. ( I90l l. we know that 

= \( V ~ 1/2 )\ ^ Iae/W e M-^(h3)m, (95) 



(AE/2Y 
m 

lo S (\\V\\-(AE/2)*) 



where 

f - , ": <96> 



Thus, we find 



IgL V exp[-di s , (< J)/a = iELfc-%,, (97) 



where we defined 

iVr = E E L < 98 > 

jeL\i iei 

dist(i,j)=r 

Note that, Ni coincides with the definition of the surface area of /, N\ = s(I). Let us 
now define the "outer boundary" of /: 

dl := {j G L\I : there exists a i G I such that dist(i, j) = 1} . (99) 
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We immediately find \dl\ < s(I) and we can restrict the sum over L\I to the set 

A r := |J {j e L\I : dist(i,i) < r - 1} , |A r | < s(i> d>r , (100) 



i.e., 



* = E E^E £i=£i*(j)I 

j£A r iel j£A r i£L j£A r 

dist(i,j)=r dist(i,j)=r 

< cr d_x Vd jf .s(/) < c 2 r %d - x s{I). (101) 



To summarize, 

log(2)(A£/2)2 



El < 4||F||c2 s m Ve-^r 2 ' 1 - 1 



r=l 

4||y||c 2 L» 1 . 2rf (e- 1 ^) 
= log(2)(Aff/2) 2 S(7) ' (1 ° 2) 

where £ii_2d is the polylogarithm of degree 1 — 2d. □ 



4 Discussion and examples 

In this section, we discuss a few special cases to exemplify the previous results. It is 
important to note that the previous results hold true in case of general lattices, beyond 
structures with a very high translational symmetry. 

Example 1 (Disordered one-dimensional system) This example is to highlight that 
besides locality, no assumptions on the coupling are required. In particular, we allow 
for random coupling, as for example in the following sense. Consider a sequence of 
one- dimensional systems on a one-dimensional chain with periodic boundary condi- 
tions, such that G^ n > = (L^ n \E^) is characterized by — [1, ...,n] and the 
adjacency matrix with entries 

E ij = + s i,j-i + &i,n8j,i + Sj >n 5i t x, i,j = 1, ...,n. (103) 

The coupling is specified by = (G^ n \ V^ n > , 1). Now let (ri, r n ) be a vector of 
realizations of random numbers taken from the interval [0, 1]. Now let (V^ n ')i t i = 3, 
and 

( V(n) )i,3 = + 1 + S i,j-l) r i + ( S i,nSj,l + Sj,n$i,l) T„, i,j = 1, ...,Tl. (104) 

From Gershgorin's theorem we then know that AE = 2A^{ 2 (V r i™') > 2. Also, 
\\V \\ < 5 for all n. Hence, Theorem 1 can be applied, and we find exponentially 
decaying correlation functions, even in this disordered system. Note that this example 
can also be generalized to d- dimensional general lattices, as one can straightforwardly 
find upper and lower bounds to the spectral values ofV using the same arguments as 
above. 
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Example 2 (Thermal states of rotating wave Hamiltonians) Consider the coupling 
C {n) = (G( n \v( n \v( n 1) on a general graph = , E^ ). Such Hamiltoni- 
ans correspond to the case of " rotating-wave Hamiltonians". For each n, the ground 
state is easy to find: it is the product state of uncoupled degrees of freedom. The co- 
variance matrix for the zero temperature case is then given by 

7 = 1 © 1 (105) 

for all n. Indeed, the correlation length is zero. For finite temperature however, we find 



7(T)=fl+2(exp(2V< n) / r )-l) J © + 2 (exp(2F (ra) / T ) - l) j. 

(106) 

Now, e.g., choose 

V {n) = t-cE {n \ (107) 

where < c < 1/2 and as in the previous example. Due to the circulant structure 
ofV M , we have AE = 2y/l - 2c, \\V\\ =l + 2c, and 

2 " e 27rfc(i-j)/n 
{XiXj) = {Pipj) = 6 itj + - ^ e 2X k (V)/T_ l - (108 > 

k— 1 

Subsequently, we merely sketch the argument leading to the statement that the corre- 
lation length is non-zero for T > 0, but it should be clear how this can be made a 
rigorous statement. For high temperature we can approximate exp(2A&(V)/T) -1« 
2Afc(V r )/T, and thus 

(xiXj) = {pifc) w 5 itj + TdV^y 1 )^. (109) 
The explicit inverse ofV^ is known A30V and given by 

t^""'" 1 '" = -w-i) ('"'"'" + """")■ ,110) 

2c 

So, z« f/ze //m/f rt — > oo, n odd, we get for i =/= j 

1 + q 2 

(iiXj) = (pipj) = — g2 Texp(-dist(i, for T > 1, (112) 

w/jere 

e = i 7T-/-T > Q- (H3) 

Hence, the thermal states have longer-ranged correlations as compared to the ground 
state. This can be explained as the first excited state is already no longer a product 
state, unlike the ground state itself. 
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Example 3 (Energy gap from exponential decay) Here we will exemplify the argu- 
ment leading to a statement on the existence of a spectral gap for a simple one- 
dimensional system. To this end consider again the one dimensional chain G^ n ' — 
(L ("),£(")) with Z>) andE^ as above and = (G("\ 1), i.e., AE^ = 
2 A^,{„ 2 (y Now suppose that in the large system limit, n — > oo, n odd, either 

(xiXj) = K cxp(-dist(i,j)/£,) (114) 



(piPj) = Kexp(-dist(i,j)/£), (115) 

with constants K > and £ > 0. We can then show that the system is gapped in 
the following way: Consider the matrix A^ 1 (cp. the previous example) the entries of 
which are given by 

That is, if All 41 holds, we have in the limit of large n 

(x l x 3 )=K^^(A~\ J , (117) 

with < q = cxp(— 1/£) < 1. Under ill 51 the same holds for (pipj). As we know 
that the inverse of A is a circulant matrix, A = 1 — cE^ n \ c = q/(l + q 2 ), we know the 
minimal and maximal eigenvalue of A are given by \\A~ 1 \\ = l + 2c, \ m in{A) = 1 — 2c 
for all n. And thus, under M14\ 

AE = 2\,„ in (V 1 / 2 ) 



\ max (v-v 2 ) 

2(1 + q 2 ) 



K(l - q 2 )\ max {A-^) 
2(1 ~q? 

Similarly, in case that ill 51 holds true, we arrive at 

t 
<1 



(118) 



AE = 2\ mbl {V 1 l 2 ) = 2K X —\\ min {A- 1 ) 

l + <r 

2K(l-q 2 ) 



(1 + 9) 2 ' 

We can hence conclude in both cases that we have that the system is gapped, 



(119) 



AE > 2(l-ex P (-l/y > q 

1 - exp(-2/£) 
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Example 4 (Algebraically decaying correlation functions) Consider a general lat- 
tice G^ of dimension d equipped with the coupling C^ n > = [G^ n \ V^ n \ 1), where 
we define the translationally invariant coupling V^ n ' via, rj > d, 

, , , „ , , f dist(i,j)~' n fori^j, 

(V^)- 1/2 =W^\ <W ln) )iJ = \ l + j:Udist{i,j)-» for, j. (121) 

Again, from Gershgorin 's theorem, we know that A mi „ (W^ ) > 1 and, together with 
the definition of the lattice dimension, Eq. 0, 



L dist{i,jY 



< l + 2j2r- r '\S r (i)\ 



r=l 
oo 

r=l 

= l + 2cC(l-d + ?7) <oo, (122) 

/or a/Z n. 77ie system is thus gapped and \\V^ n '\\ < 1, i.e., all assumptions of Theorem 
1 are met except that y( n ) is not of finite range. Hence we cannot conclude that the 
correlation functions are exponentially decaying. Indeed, we find that the position 
correlation function for i =/= j is given by 

^ = msihv' (123) 

so the correlations in this non-local system are algebraically decaying, despite being 
gapped. 



5 Outlook 

In this paper, we have rigorously revisited the question of algebraically and exponen- 
tially decaying correlation functions in gapped harmonic systems on general lattices. 
Both for the ground states, as well as for Gibbs states, we considered the implications 
of the gap and the coupling strength to the correlation length. For systems only coupled 
in position, we showed that an energy gap can be deduced from algebraically decaying 
correlations. We also found an equivalence between the existence of a gap and expo- 
nentially decaying correlation functions for local Hamiltonians. For zero temperature, 
no assumptions have been made on the underlying lattice. This harmonic quasi-free 
case is a particularly transparent instance of a physical system where this physically 
plausible connection can be explored, explicitly making use of properties of matrix 
functions of banded matrices, instead of exploiting Lieb-Robinson type results. As 
such, we cover cases not included in the class of systems considered in Ref. |4). 
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We showed that these statements have immediate implications on the scaling of en- 
tanglement in ground and thermal states of quantum many-body systems: under these 
conditions, the geometric entropy, so the von-Neumann entropy of the reduction, is 
bounded from above by a quantity that is linear in the boundary area of the distin- 
guished region. Equivalently, this quantifies the degree of entanglement of the region 
with respect to the rest of the lattice. This result on a connection between the surface 
area and the degree of entanglement can hence also be identified on general graphs, 
hence establishing an rigorous entanglement-area relationship in general finite-ranged 
gapped harmonic lattice systems. 
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